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Abstract. - By large-scale Monte Carlo simulations of semiflexible polymers in d = 2 dimensions 
the applicability of the Kratky-Porod model is tested. This model is widely used as "standard 
model" for describing conformations and force versus extension curves of stiff polymers. It is 
shown that semiflexible polymers in d = 2 show a crossover from hard rods to self-avoiding walks, 
the intermediate Gaussian regime (implied by the Kratky-Porod model) is completely absent. 
Hence the latter can also describe force versus extension curves only if the contour length is only 
a few times larger than the persistence length. Consequences for experiments on biopolymers at 
interfaces are briefly discussed. 
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Characterizing the flexibility or stiffness of polymer 
chains is of basic importance for describing their struc- 
ture and dynamics, and hence relevant for understand- 
ing the functions of biopolymers, as well as the applica- 
tion properties of synthetic polymers ]T]-j3]. Moderately 
stiff ("semiflexible") macromolecules behave like rods on 
small scales, and one captures this behavior by the con- 
cept of the so-called "persistence length" ip. For larger 
length scales, entropic flexibility prevails and random coil- 
like structures occur. Important examples for such stiff 
biopolymers are DNA, some proteins, actin, neurofila- 
ments, but also mesoscopic objects such as viruses [5- 7 . 
The experimental study of such biopolymers and the inter- 
pretation of these observations by models is a very active 
topic of research {e.g. [5HT7]}. In particular, the confor- 
mation of these biopolymers can be directly visualized by 
electron microscopy (EM) or scanning force microscopy 
(SFM) techniques when such polymers are adsorbed on 
substrates f8HIO l [T2 j [T4 j [T7]: by atomic force microscopy 
(AFM) also force versus extension curves can be mea- 
sured [TT1[T3]. The same methods also work for synthetic 
polymers such as molecular brushes [18] , where stiffness is 
controlled by the length of side chains [19]. 



The standard theoretical model, that is almost exclu- 
sively used {e.g. [20H3T] } to interpret these experiments 
is the simple "wormlike chain (WLC) model" [Mll33]. Its 
Hamiltonian is, in the continuum limit. 
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Here the curve f[t) describes the linear macromolecule, 
t is a coordinate along its contour which has the length 
L. We choose units such that ksT = 1, and the bending 
stiffness k then is k = in d = 2 dimensions. In this 

paper we shall focus on the case of chains confined to two- 
dimensional geometry, since this case is relevant for the 
EM and SFM imaging techniques, and also the subject of 
numerous theoretical studies (e.g. [25 1I28H30] ). However, 
the applicability of eq. ([Ij in principle is questionable, 
since it neglects excluded volume between the repeat units 
of the chain completely. Thus, eq. (|T]) yields the end-to- 
end distance of the polymer chains as 



= 2^pL|l-i[l-exp(-n)]| 



L/lp , (2) 
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and hence for n 3> 1 the chain behaves like a Gaussian 
coil ((i?^) = 2lpL) while for n < 1 it is essentially a rigid 
rod of length L. The bond-autocorrelation function shows 
then a simple exponential decay, 

g{t) = (a, • a^+s) = exp(-i/^p) , t ^ sib , (3) 



where we now consider a chain where bonds of length 
£{, connect repeat units at sites f^, = r^+i — fi, \ai\ = £{,; 
so L = N^lb- Finally, if one considers the effect of a force 
/ acting on one chain end (the other being fixed at the 
origin), by adding a term - fX to the Hamiltonian {X 
being the x-component of the end-to-end distance), one 
obtains from eq. ^ the force vs. distance relation to a 
very good approximation, in d = 2 [25 



(i?2)l/2 0C^y*i3/4 ^^^^^ 

and g{t) oc t"^/^, for t > £p, rather than eq. The 
latter result is consistent with the scaling prediction [36] 
g{t) oc with (3 — 2(1 — where the Flory exponent 
1/ = 3/4 in d = 2, as written already in eq. ([5|). 

There has been evidence for the scaling (i?^)^/^ oc L^^'^ 
for not so stiff polymers such as single stranded DNA in 
d = 2 dimensions, see e.g. PllOlllGj . but it has been widely 
believed that for very stiff polymers excluded volume in- 
teractions (that cause the nontrivial exponent v — 3/A 
rather than the Gaussian result i' — 1/2 which follows 
from eq. ([2])) can be neglected, except for extremely long 
chains. We will show, however, that excluded volume ef- 
fects set in strongly already for L w 5£p, invalidating the 
straightforward use of eqs. ©-(jl]) for many cases of inter- 
est. 

We carried out Monte Carlo simulations of self-avoiding 
walks (SAWs) on the square lattice, applying an energy et 
if the orientation of bond vector differs (by ±7r/2) from 
that of di-i, and using the pruned-enriched Rosenbluth 
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Fig. 1: Semi- log plots of the bond-correlation g{t) vs. the con- 
tour length t, for the ranges (for a definition of the parameters 
see eq. ^ 0.1 < qb < 1.0 (a) and 0.005 < qb < 0.05 (b). The 
data are taken for L — 25600 and 6=1, averaging over the 
site i in eq. Straight lines indicate fits of the initial decay 
of g{t) to eq. ([Sj. The resulting values of £p are quoted in the 
figure, (c) Log-log plot of g(t) vs. t, for qt = 0.005 to 1.0 (from 
above to below). The straight line shows a fit of the data for 
qb ~ I and f > 10 to the power law g(t) oc t""'^. 



/^. = i[6^-i+(i-^n • (4) 

Since various experimental data have been described by (b) 
eqs. ©-dl]) with some success adjusting parameters such 
as £p and L, it is widely believed that the basic Kratky- 
Porod model, eq. ([T]), describes semiflexible chains accu- 
rately, and a large body of work is concerned with various 
refinements of this model {see e.g. [5BH30]}- 

However, in the present Letter we show that in fact in 
d = 2 the validity of the Kratky-Porod model in the good 
solvent regime is very restricted, it always holds only up 
to contour lengths L of a few times £p, irrespective how 
large the persistence length £p is. In particular, a regime 
of L where Gaussian statistics holds, (i?^) — 2£pL, in d — 
2 is completely absent, unlike the case of d = 3, where 
for very large £p a double crossover (rods — > Gaussian 
coils — > non-Gaussian swollen coils) is established both (^q^ 
experimentally [33] and theoretically [5S]. Also eq. (|3]) 
breaks down for L ^ £p, irrespective how large £p is. In 
d — 2, we will show that 
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method [3S1[371[35] . The partition function of SAWs with 
Nt, steps and A^bcnd local bends is 

config 

where qt — e:sq){—eb/kBT), b ~ cxp{f /ksT) and X is the 
x-component of the end-to-end distance (assuming that 
the force / acts in the -|-a;-direction). In experiments 
where a force is applied to an end of a strongly adsorbed 
chain, that takes essentially two-dimensional conforma- 
tions, it is possible to direct this force either perpendicular 
or parallel to the surface; only the latter case is considered 
here. Note qb — ^ for flexible chains (standard SAWs) and 
6 = 1 in the absence of the force /. We generated data for 
C(7V6, A^bend,^) for 0.005 <qb< 1.0 and Nb < 25600. 

Fig-IDshows the bond-orientational correlations (for the 
case / = 0). For rather flexible chains, qb = 0.4, there are 
at best a few values i = 1, 2, 3 compatible with an expo- 
nential decay (we use ib = 1 here and in the following) . For 
small qb, eq. ([3]) has a more extended range of applicability, 
and £p strongly increases when qb decreases, £p ~ 0.61/gb. 
But the asymptotic decay always is the expected power 
law (fig. HJc)). As has been emphasized recently |39j, in 
the presence of excluded volume "the" persistence length 
is a somewhat ill-defined concept; for the present model, 
ip henceforth is defined from the initial slope of the curves 
\ng{t) vs. t as t 0. 

Fig. m presents a test of eq. While eq. ^ trivially 
works for L < £p (the rod- like regime), significant devia- 
tions become visible for L > 5ip, irrespective of how large 
£p is, as the scaling plot (fig.lljb)) shows. In contrast to oc- 
casional claims in the literature |12] , a regime of Gaussian- 
like coils is completely absent in d = 2. This result can 
be rationalized by the proper adaptation of Flory-type ar- 
guments |40] to d = 2. The free energy of a stiff chain is 
taken as the sum of an elastic energy [B? /£pL) and the 
enthalpy due to repulsions, proportional to the 2nd virial 
coefficient {v2 = £p [H] ; prefactors of order unity are sup- 
pressed throughout) 
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Fig. 2: Log-log plot of (i?^)/(2L) versus L = Ndb (a) and 
log-log plot of {B?)/{2£pL) versus L/£p (b), for h—1 and sev- 
eral choices of qb, as indicated. Full curves show the WLC 
prediction, eq. ((2}, using Ip (highlighted by arrows in (a)) 
from fig. [ira)(b) as an input. Straight lines in (b) indicate 
the power laws in the rod regime {{Ft?) ~ L?) and the SAW 
regime {eq. (O}, respectively. 



A^^ = Ry{£pL) + v2R'[{L/£p)/R^f 



(7) 



In d = 2, the "volume" of a chain of radius R scales like R^ , 
and the density of the n = Ll£p sub units is n/R^ in this 
volume. Minimizing IS.F with respect to R yields eq. ([5]). 
The minimum length L where eq. ((S]) holds is found when 
the enthalpic term in eq. Q is unity for R? = £pL, i.e. for 
L* = £p/v2 = £p, and there the rod-like regime starts: this 
argument shows that we should expect a single crossover 
from rods to SAWs, as seen in fig. [Sfb), unlike the d — 3 
case [351140]. 

How then can we understand the apparent success (sug- 
gested in the literature) of the Porod-Kratky model to an- 
alyze force-extension curves in 2d? In fig. [3]we show some 
of our results on force vs. extension curves in d = 2 and 
compare our data to the theoretical prediction based on 



the WLC model, eq. (jH) . Here the persistence length esti- 
mates quoted in Figs. [ira)(b) were used, so we can com- 
pare our simulation results that are based on eq. © to 
the prediction, eq. ([4]), without adjusting any parameter 
whatsoever. One can see that the latter equation works 
only approximately (fig. EI a)) for very stiff chains in a very 
restrictive range of contour lengths, where we can deduce 
from a detailed inspection of the data that 6 < L/£p < 10 
must be fulfilled: if L/£p is too small, the chain behaves 
as a flexible rod, which can be oriented by a force but 
not stretched; if L/£p is too large, excluded volume ef- 
fects invalidate eq. @, similarly as eq. ([2]) fails then. For 
qb — 0.4, the chains have hardly any rod-like regime as 
Fig.IHa) reveals, £p is less than three lattice spacings, and 
so large deviation from eq. (|4]) are no surprise, of course. 
For small qb, where for the chosen values of L = 200 in 
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Fig. 3: Log-log plot of {X)/L vs. f£p for L = 200. Eq. ^ is 
shown by the full curve for comparison (a) . Rescaled force f£p 
plotted against {X)/L for L = 1600 (b) and L = 25600 (c). 
Various values of qt are shown as indicated. 



Fig.[3I^a) L is only a few times larger than ip (recall £p ~ 62 
for qb — 0.01, Fig. [Hb)), the deviations of the data from 
eq. dl]) go into the opposite direction {{X)/L for f£p < 1 
is smaller than predicted by eq. (|4]), while {X)/L is larger 
than predicted if £p is small). This finding implies that 



for L ~ 200 and intermediate values of ip, the observed 
variation of {X)/L with f£p is close to the predicted one, 
for the intermediate range of L/£p quoted above, but this 
agreement is somewhat accidental. 

Note also that a sensible test of the Kratky-Porod model 
(which is a continuum model) by our discrete lattice model 
is only possible for forces such that f£p < 1, since im- 
portant deviations between discrete chain models and the 
Kratky-Porod model occur [551 when the so-called deflec- 
tion length A oc (/^p)~^^^ of worm- like chains becomes 
smaller than the bond length £i,. Thus our data do not 
converge to eq. ^ even for large f£p, although for very 
strongly stretched chains {{X)/L) close to unity) excluded 
volume effects must become irrelevant. 

If L is very large, such a crossing of the simulated curves 
for (X) /L as function of f£p with cq. ^ when £p is var- 
ied does no longer occur (Fig. [31Ib)(c)). The simulation 
results for (X) / L are now always significantly larger than 
the prediction, eq. ([4]), particularly for small values of fip. 
This huge discrepancy for small values of f£p can be un- 
derstood readily in terms of a linear response argument: 
actually, eq. (jlj is found from adding a term —fX to the 
Hamiltonian, eq. ([T]). Therefore it is straightforward to 
derive, in the limit / — 0, the linear response relation 



(8) 



Since (X^)/=o — (i?^)/2, where according to the Kratky- 
Porod model {eq. (I2|)} for L 3> £p we have simply (i?^) = 
2£pL, we conclude that {X) = {X'^)f = £pfL (in agree- 
ment with the Taylor expansion of eq. (|4]) to first order in 
{X)/L, as it must be, of course). However, in d = 2 for 
vanishing force and L ^ £p the relation (X^) — £pL must 
be replaced by (X^) cx £l^'^L^/^, as is readily seen from 
eq. ([5]). Therefore we predict for the linear response regime 
a very different scaling for the force-extension behavior, 
namely 

{X)/L^£l^'L'/\f. (9) 

This relation is tested in Fig. |4l A wide range of choices 
of contour lengths L = Ni,£i, and several choices of qb 
and hence £p (the relation between qb and £p is quoted 
in Fig. [IJa)(b)) are included. An interesting issue also is 
the regime of relative extensions over which linear response 
holds: while eq. (j4|) implies a linear response regime apply- 
ing almost up to {X)/L « 0.3, irrespective of £p, we sug- 
gest that the linear response breaks down if {Xy « {X"^), 
i.e. for {X)/L cx {£plLfl^ ^ as £plL 0. In the 
nonlinear regime, fig. I31[b)(c) suggests that {X) / L can be 
described by some universal function oi £pf, that does not 
depend on £p: this is the universality of c? = 2 SAWS, not 
the Kratky-Porod model. 

Of course, the scaling (X) oc L^^^f for small / is con- 
sistent with the scaling behavior proposed by Pincus [42] 
for stretched flexible polymers in the presence of excluded 
volume 

(X) = RoF{Ro/^p) (10) 
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Fig. 4: Log-log plot of f{LlpY^'^ vs. {X)/L, including several 
values of gj, as indicated, and data for Nt = 400, 1600, 6400, 
and 25600 (from bottom to top at the right side of the diagram, 
respectively). Straight dotted line (black) indicate the linear 
response behavior, and straight solid line (red) indicate the 
non-linear behavior, i.e. / oc (X)^ (see text). 



where Rq is the radius of chain in the absence of a stretch- 
ing force, F{Ro/^p) is a scaling function, and oc 1// 
is the size of "Pincus blobs", and hence in the linear 
response regime (X) oc Rq/, i.e. eq. ^ results. The 
condition that {X)/L is of order unity then leads to [52] 
{X} oc f^^"^^ — /^/^ in d = 2 dimensions, i.e, a strongly 
non-linear relation between / and {X). This power law in 
fact is compatible with the data in Fig. 2] for large enough 
{X)/L. 

In conclusion, we have shown that in d = 2 dimensions 
the Kratky-Porod model, that is ubiquitously used to ana- 
lyze the internal end-to-end distances of biopolymers such 
as DNA {e.g. [HIIIT]} or of synthetic polymers such as 
the bottle brushes {e.g. [18]} and to analyze force versus 
extension curves {e.g. [TTlfTS] } has a very limited validity: 
it trivially describes the rod-like regime, L < ip, but this 
regime is not useful in the context of such measurements, 
which are devoted to understanding the dependence of the 
persistence length on various parameters (such as partic- 
ular amino acid sequences in DNA, or side chain length in 
bottle brushes, etc.). In d = 2, a regime where Gaussian 
statistics (requiring L ^ £p) holds is completely absent. 

Our findings imply that conformations of semiflexible 
polymers in d = 2 (equilibrated surface adsorbed case) 
depend on their relative length L/£p very differently from 
the case d — 3 (dilute bulk solution) . Thus there is no di- 
rect way to infer properties (such as ip) in the bulk from 
measurements on surface adsorbed chains: there is no sim- 
ple relation between the effective persistence lengths either 
{in our model £p oc l/(4q{,) for d — 5 but £p = 0.61/gt, in 
d = 2 for gf, }. 

Going beyond the strictly 2-dimensional case, exploring 



the crossover to weak adsorption (chains with dangling 
non-adsorbed "tails" and "loops" in addition to adsorbed 
"trains") will be intriguing. Also, the effects of excluded 
volume on force versus extension curves when strongly 
adsorbed chains are pulled off a surface in the direction 
normal to the surface by an AFM tip need to be studied 
carefully. Thus, much further work is needed for a bet- 
ter modeling of biopolymers and other stiff polymers at 
interfaces, and on the interpretation of the corresponding 
experiments. 

The effects studied in our work should also be relevant 
when one studies semiflexible chains confined to the sur- 
face of a sphere or its interior [43] , a problem believed to 
be of great biological relevance. 
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